Magnetic Ftelds across the Hertzsprung-Russell Diagram 

ASP Conference Series, Vol. tt x 10 7 , 2001 

G. Mathys, S.K. Solanki, & D. T. Wickramasinghe 

Further Analysis of Stellar Magnetic Cycle Periods 

Steven Saar 

Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, 
Cambridge, MA 02138, USA 

Axel Brandenburg 

NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

Abstract. We further investigate relationships between activity cycle 
periods in cool stars and rotation to include new cycle data, and explore 
different parameterizations of the problem. We find that relations be- 
tween cycle and rotational frequencies (w C yc vs. 0) and between their 
ratio and the inverse Rossby number {uj cyc /VL vs. Ro _1 ) show many sim- 
ilarities, including three branches and similar rms scatter. We briefly 
discuss some implications for dynamo models. 



1. Introduction 

Several recent studies (Ossendrijver 1998; Tobias 1998; Brandenburg et al. 1998; 
Saar & Brandenburg 1999 [=SB]; Lanza & Rodono 1999) have revisited rela- 
tionships between stellar magnetic cycles and other stellar properties, taking 
advantage of the increased quality and quality of the cycle data available (e.g., 
Baliunas et al. 1995). SB studied relationships between non-dimensional quan- 
tities such as cycle-to-rotational frequency ratio tu cyc /VL, the normalized Ca II 
HK emission flux -R^ K , and the inverse Rossby number Ro _1 = 2t c Q (where 
r c is the convective turnover time). They found evidence for three power-law 
"branches" upon which stars tended to cluster. Here we expand on this work. 
We add new cycle data, and investigate how the new data affect various parame- 
terizations, both dimensional and non-dimensional, of the stellar cycles, focusing 
on relations between uj cyc and rotation. 

2. Data and Analysis 

We combine cycle and stellar data gathered by SB with more recent measure- 
ments of plage (e.g., Hatzes et al. 2000) and spot cycles (e.g., Olah et al. 2000). 
Cyclic changes in P ro t in some close binaries have been linked with magnetic 
cycle modulation (via changes in mean magnetic pressure) of stellar quadrupole 
moments (Lanza et al. 1998). These cycles based on variations in P rot (Lanza & 
Rodono 1999) are also tentatively included. We follow the strategy of SB, using 
theoretical r c (Gunn et al. 1998) and weighting the P cyc by a "quality factor" 
w (0.5 < w < 4) depending on the strength of the periodogram signal or clarity 
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of the cycle modulation. We set w = 1 for the P rot -change cycles. Stars are 
assigned to branches (where appropriate) by eye to minimize fit rms. Evolved 
stars were not included in the fits due to less well determined r c . Results for 
different classes of stars are shown in Fig. 1 LEFT (using dimensionless w cyc /ft 
and Ro- 1 ) and Fig. 1 RIGHT (using uj cyc =2vr/P cyc and ft = 2vr/P ro t). 

3. Results and Discussion 

Our results can be summarized as follows: 

(1) Three branches - denoted I (inactive), A (active), and S (super-active) 
- appear in both the Ro _1 and ft parameterizations (Fig. 1). For the Ro _1 fit, 
the power law exponents are 5i « —0.3 (with a fit dispersion a fit = 0.095 dex), 
5 A ~ -0.15 (a fit = 0.18), and 8 S ~ 0.4 (a fit = 0.26 dex); for the ft fit, 5 £ « 1.15 
{a ^ = 0.093 dex), 5 A ~ 0.8 (a fit = 0.17 dex), and 5 S « 0.4 (a fit = 0.24 dex). 
Thus the rms scatter is similar for the two parameterizations. 

(2) 

(2) Secondary cycle periods (-Fy t ) seen m some stars often lie on one of 
the branches (though this is more rare in S branch stars). The solar Gleissberg 
"cycle" (~ 100 years) appears to lie on the S branch. The preferred branch of 
the primary P cyc (with the strongest periodogram signal) may be mass and ft de- 
pendent. Multiple P cyc may reflect multiple dynamo modes in an aft framework 
(Knobloch, Rosner & Weiss 1981), or different dynamos existing in separate lat- 
itude zones (note the dual, separately evolving activity patterns in the double 
P cyc star (3 Comae; Donahue & Baliunas 1992). In the Babcock-Leighton sce- 

(2) 

nario, P rot may be excited by stochastic variations in the poloidal source term 
(Charbonneau & Dikpati 2000). 

(3) A single power law can be fit to the data (e.g., uj cyc c(. ft~°' 09 , SB; see also 
Baliunas et al. 1996) but only at the expense of a considerably higher dispersion 
about the fit {(J fa = 0.33 dex), and loss of an explanation for the secondary cycle 
periods (since they no longer reside on another dynamo "branch"). 

(4) Evolved stars typically lie near branches, though show more scatter 
than the dwarfs. Since the increased scatter is seen in both parameterizations, 
it is unlikely to be due to less precise r c in evolved stars (indeed, arguably the 
scatter in evolved stars is reduced using Ro _1 ). The P cyc based on P rot variation 
(Lanza k, Rodono 1999) also follow the general trends. The branches are better 
separated using Ro _1 . On the other hand, the ft plot is simpler, lacking the 
"transitional" regime between the A and S branches seen in the Ro _1 diagrams. 
Contact binaries (gray o; bottom panels) are poorly fit in both schemes (worse 
if Ro" 1 is used); their dynamos may be altered by turbulent energy transfer 
toward the secondary (Hazlehurst 1985) which is independent of rotation. 

(5) The branches may merge for small Ro^ 1 or ft (though at values which 
might not be reached by actual stars). Curiously, the ratio of the power law 
exponents for the ft fits are 5/ : 5a '■ 3s ~ 3 : 2 : 1. 

(6) Since ft and Ro _1 decrease in time on the main-sequence, the relations 
between uj cyc and rotation map out dynamo evolution with time. The overlap- 

(2) 

ping branches and P rot suggest that uj cyc evolves in time in a complex, sometimes 
multi-valued fashion. The panels of Figure 1 LEFT show an approximate age 
calibration along the top axes. 
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A Babcock-Leighton type model predicts lu cjc (x u^ 9 for solar-like dwarfs 
(where u m is the meridional flow velocity; Dikpati k Charbonneau 1999). If u m 
increases approximately linearly with Q in slower rotators (e.g., Brummell et al. 
1998), the predicted u) C y C matches the I and A branches reasonably well (see also 
Charbonneau k Saar, this volume). Mean- field models with sufficiently strong 
£1 dependence for the differential rotation (e.g., Donahue et al. 1996) and the a 
effect (e.g., Brandenburg k Schmitt 1998) can also match the observed branches 
(SB; Charbonneau k Saar, this volume). We are studying a variety of dynamo 
models to better understand the implications of the cycle - rotation relations 
seen here. 
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Figure 1. LEFT top: w cyc /f2 vs. Ro 1 for single dwarfs; symbols indicate 
the sun (0), F (A), G~fo), and K (□) stars (filled if log i?^ K > -4.75; size 
oc ^/w, the P C yc "reliability"). Dotted vertical lines connect two P cyc (a x 
marks frot)> or ^cyc with a long-term trend (i.e., a possible P^ot > 25 yr; 
arrow symbol). Weighted least square fits (cj cyc /Sl oc Ro 5 ) for the active (A) 
and inactive (I) branches are shown (solid); Sj = —0.32 and 5a = —0.16. 
LEFT middle: same, including binaries (BY Dra, CV secondaries; M stars = 
o) and RS CVns (+; not included in the fits). A "superactive" (S) branch 
appears, with 5s = +0.43 (P rot -Pcyc lines shown only for new stars). A 
transitional regime between the A and S branches is indicated (dash-dot). 
LEFT bottom: same, including cycles based on P rot variation in RS CVns 
(new +), CV secondaries (open o), Algols (*), and contact binaries (gray o). 
RIGHT top: uo cyc vs. f2 for single dwarfs. Fits (u} cyc oc Q s ) for the A and 
I branches (solid) yield 5i = 1.15 and 5a = 0.80. RIGHT middle: same, 
including binaries (like LEFT middle). The new S branch shows 5s — 0.38. 
RIGHT bottom: same, including P rot variation cycles (like LEFT bottom). 



